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O ■ We study the semi-holographic idea in context of decaying dark components. The energy 

i-C ' flow between dark energy and the compensating dark matter is thermodynamicaUy general- 

^ . ized to involve a particle number variable dark component with non-zero chemical potential. 

^ ' It's found that, unlike the original semi- holographic model, no cosmological constant is 

' needed for a dynamical evolution of the universe. A transient phantom phase appears while 

a non-trivial dark energy-dark matter scaling solution keeps at late time, which evades the 

^ . big-rip and helps to resolve the coincidence problem. For reasonable parameters, the deceler- 

O 

' ation parameter is well consistent with current observations. The original semi-holographic 

m 



model is extended and it also suggests that the concordance model may be reconstructed 
^SJ ! from the semi-holographic idea. 

(N 



PACS numbers: 98.80.-k 95.36.-^x ll.lO.Lm 

I. INTRODUCTION 

The existence of dark energy is one of the most significant cosmological discoveries over the 
last century [1]. Various models of dark energy have been proposed [2], such as a small positive 
cosmological constant, quintessence, k-essence, phantom, quintom, etc., see Ref.j^ for a recent 
review. Stimulated by the holographic principle, it's conjectured that dark energy problem may be 



a problem of quantum gravity, although its characteristic energy scale is very low 



According to 
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the holographic principle^], the true laws of mside any surface are actually a description of how its 
image evolves on that surface. In the case of black hole, the event horizon is a proper holography. 
When applying the holographic principle to the universe, the event horizon is probably not a good 
candidate of holography screen at least for two reasons: the identification of the cosmological 
event horizon is causally problematic in that it depends on the future of the cosmic evolution; a 
decelerating universe would have no holographic description, since it has no event horizon at all! 
With reference to the seminal work of Jacobsony], several arguments have been put forward that 



the apparent horizon shou 
and Hawking temperature 



dd be a causal horizon and is associated with the gravitational entropy 
7 1 Q], and hence the right holography of the universe. Apart from this 
key observation, other thermodynamical properties of the dark energy and dark matter have yet 
to be found in detailsyj. Recently, the semi-holographic dark energy model[l^ has been proposed 
to examine how the universe could evolve if a dark component of the universe strictly obeys the 
holographic principle. It was found that, based on the first law of thermodynamics, the existence of 
the other dark component (dark matter) is compulsory, as a compensation of dark energy. In that 
model, with a non-zero cosmological constant the standard cosmological kinematics is recovered and 
the cosmological coincidence problem is also alleviated with the existence of a stable Einstein-de 
Sitter scaling solution. 

For all that, the properties of the semi-holographic model should be clarified for further estima- 
tions from three points of view. First, only massless particles (photons) exchange are considered. 
Or we treat dark matter and dark energy as two closed systems. Surely, by nature they can be 
open systems, and thus the particle numbers of the two systems can be variables. Second, the EOS 
of dark energy is constrained to be LOde < — 1. It certainly is in accord with the latest cosmological 
observations and phenomenologically reasonable. However, as a thermodynamics motivated dark 
energy model, the construction would be founded on a more solid basis of thermodynamics and 
statistical physics, and exactly in this sense, the well-known theoretical difficulties does matter. As 
a matter of fact, the case of phantom energy with w < — 1 is ruled out because the total entropy 
of the dark component is negative[lll]. The equation of state (hereafter EOS) was restricted to the 
interval — 1 < uj < —1/2 and a fermionic nature to the dark energy particles is favored. Thermody- 
namics arguments [is] in favor of the phantom hypothesis had to resort to an unusual assumption 
that the temperature of a phantomlike fluid is always negative in order to keep its entropy positive 
definite (as statistically required) or by arguing that the scalar field representations of a phantom 
field has a negative kinetic term which quantifies the translational kinetic energy of the associated 
fluid system. Third, as is known, the possibility of a coupling between dark matter and dark energy 
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is often considered with three types in the hterature, and those interaction forms are always intro- 
duced phenomenologically or through the low energy effective action [l^: dark matter decaying into 



dark energy 



12], dark energy decaying into dark matter [l3|] and interaction in both directions 14 1. 



It's worth noting that, due to the non-adiabaticity of the holographic evolution of dark energy, 
the semi-holographic model obviously introduced a dark matter/dark energy non-gravitational in- 
teraction (A hypothetical dark energy prop erty sometimes induces a non-gravitational coupling 



between dark matter and dark energy. See 



16l | for another interacting model of this kind.), and 



the interaction belongs to the last category: dark matter dominated in the past and transferred 
energy to dark energy; the latter begins to dominate at present, and in the near future, the transfer 
would be reversed and dark energy decays to dark matter with a stable scaling solution reached. 
The theoretical obstacle is that, provided that the chemical potential of both components vanish, 
it is the decay of dark energy into dark matter that was favored by the thermodynamical point of 
view and the second law. 

Owing to a series of systematic work[ll| {l?! [l^ [lol j2o| . thermodynamics and statistical argu- 
ments on dark components and the cosmic evolution have come to several important conclusions 
which may kick off those dilemma. By considering the existence of a non-zero chemical potential, 
reanalysis of the thermodynamics and statistical properties of the dark energy scenario supports 
that the temperature of dark energy fluids must be always positive definite and that a nega- 
tive chemical potential will recover the phantom scenario without the need to appeal to negative 
temperatures [19]. In addition, a bosonic nature of the dark energy component becomes possible. 
If the chemical potential of at least one of the cosmic fluids is not zero, the decay can occur from 
the dark matter to dark energy, without violation of the second law 20|]. From this point of view. 



the decay process between dark components with non-zero chemical potential is not merely a the- 
oretical tool, but should be seriously considered as a prerequisite to turn the phantom dark energy 
and the decay in both directions to be physically real hypotheses. 

More importantly, due to the insufficiency in discussing the particle number variable transferring 
process between dark components, the original semi-holographic model should not be regarded as 
a complete thermodynamical examination of the semi-holographic idea in cosmology. A general 
framework dealing with the decay of the dark component particles is still in lack and therefore 
well motivated. Particularly speaking, the first law of thermodynamics would be applied in the 
context of grand canonical ensemble rather than canonical ensemble describing the original semi- 
holographic model. 

This paper is organized as follows: In the next section we will study the semi-holographic 
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model in grand canonical ensemble with the dark energy chemical potential explicitly given. The 
dynamical analysis is left in section III. We find that there exists a stable dark matter-dark energy 
scaling solution at late time, which is helpful to address the coincidence problem. After a close 
analysis and comparison between model behavior and the observational facts, we present our 
conclusion and some discussions in section IV. 



II. THE MODEL 



There is decisive evidence that our observable universe evolves adiabatically after inflation in 
a comoving volume, that is, there is no energy-momentum flow between different patches of the 
observable universe so that the universe keeps homogeneous and isotropic after inflation. That is 
the reason why we can use an FRW geometry to describe the evolution of the universe. In an 
adiabatically evolving universe, the first law of thermodynamics equals the continuity equation. In 
a comoving volume the first law reads, 

dU = TdS - pdV + ildN, (1) 

where U = ^kP<^'^ is the energy in this volume, T denotes temperature, S represents the entropy of 
this volume, V stands for the physical volume V = ^kO? , jl is the chemical potential of the energy 
component with particle number N . Here, is a factor related to the spatial curvature. For 
spatially flat case VLq = |7r, in this paper we only consider the spatially flat model, p is the energy 
density and a denotes the scale factor. The last term in the above equation indicates that the grand 
canonical ensemble instead of the canonical ensemble for the original semi-holographic model is 
considered. As a consequence, the framework allows decaying of the particles in consideration. 

The semi-holographic model concerns the possibility of a non-adiabatical dark energy, where 
the term TdS does not equal zero. Based on the investigations in [3, ^, 
volume the entropy becomes. 



the entropy in a comoving 



be- ^^jp^ = 8tt ^1 Ha . (2) 

where H is the Hubble parameter, /i denotes the reduced Planck mass. The entropy has been 
reasonably assumed to be homogeneous in the observable universe. As our observable universe 
evolves adiabatically after reheating, the varying entropy of dark energy in a comoving volume 
should be compensated by an entropy change of the other component to keep the total entropy 
constant in a comoving volume. 



With the above supposition and conventions the entropy of the dark energy satisfies ^ 



Sde = STT^H^Ha^. 



(3) 



Correspondingly, the entropy of dark matter in this comoving volume should be 



'dm 




(4) 



where C is a constant, representing the total entropy of the comoving volume. 



The Friedmann equation will determine the evolution of our universe 



H'^ = ir^iPdm + Pde + A), 



(5) 



where H is the Hubble parameter, denotes the density of non-baryon dark matter, p^e denotes 
the density of dark energy, and A is the cosmological constant (or vacuum energy). Partly because 
the partition of baryon matter is very small and does little work in the late time universe, and 
partly because the non-gravitational coupling between baryonic matter and dark energy is highly 
constrained, we just omit the baryon matter component in the discussion. The appearance of the 
cosmological constant in the expression is just for a general discussion, and as can be seen below, 
it is not a necessary component in the generalized semi-holographic model. As a comparison, the 
cosmological constant has to be present in the original model, or else the ratio of holographic dark 
energy and dark matter density would always remain constant which is obviously not realistic. 



with the efficient a a positive constant. As was shown by Pereira, if the coefficient d > 0, a negative 
chemical potential will make the phantom fluid thermodynamically consistent and a decay of dark 
matter into dark energy possible. This assumption may help us highly broaden the exploration of 
the cosmological parameter region motivated by thermodynamics. 



Below we will assume the chemical potential of dark energy to be of the form in Ref.[20] 




jl = —aT 



(6) 



III. DYNAMICAL ANALYSIS 



To investigate the evolution in a more detailed way, we take a dynamical analysis of the universe. 
The holographic principle requires that the temperature j3] 
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By using ([7|), ([5]), and using the particle number A'" to be proportional to the energy density, 
the first law of thermodynamics ([I]) becomes the evolution equation of dark energy, 

^ / Pdm 

3^'^ TT&Y ' 

where a prime denotes the derivative with respect to In a, indicates the EOS of dark matter. 



Wde represents the EOS of dark energy and Y = \J pde + Pdm + A. Similarly, we derive the evolution 
equation of dark matter, 



2 , _ pde[-l + (1 - 2aY)wde] - Pdm[3 + 2aY + (1 + 2aY)wdm] - 2A 

Clearly, the above equations will degenerate to the equations (8) and (9) of the paper [l^ when 
d = 0. For convenience we introduce two new dimensionless functions to represent the densities. 



A Pdm f-ir\\ 

^ = :J72 772 ' (10) 



(11) 



A Pde 

and a dimensionless cosmological constant 

where Hq denotes the present Hubble parameter. Then the equation set (l9|), ([8]) becomes 
2^, ^ -u[3 + 2V3ay + (1 + 2V3aywdm)] + ^[-1 + (1 - 2V3ay)wde] - 2A 

2 , _ u{l - Wdm) - vj'^ + 2^/3ay + 3wde) + 2X , , 

S"" ~ 1 + V3ay ^ ^ 

respectively, where y = \/u + v + X and a = dpHo are both dimensionless. We note that the time 
variable does not appear in the dynamical system (|13p and (|14p because time has been completely 
replaced by scale factor. 

Before presenting the numerical examples for special parameters we study the analytical prop- 
erty of this system. And we will below take A = for simplicity. The critical points of the 
dynamical system and (fn|) are given by 

< = v', = 0, (15) 

which yields, 

Uc = 0, Vc = (16) 
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and 

(1 + Wde)[l + Wde{2 + Wdm)? [1 + W^de(2 + Wdm)? 



(17) 



a'^iWde-Wdm){'i-+Wdmy ' a2(ti;rfe - + Wdm) ' 

The trivial solution corresponds to a future infinitely diluted universe and is less interesting 
for the coincidence problem. On the other hand, the non-trivial fixed point satisfies the below 
equality: 

Vc 1 + Wdm ' 

which is shared by the seminal model without chemical potential involved. So, finally the universe 
enters a de Sitter phase, and the ratio of dark matter over dark energy is independent of the chemical 
potential. Furthermore, there are two reasonable cases for the non-trivial scaling solution: case I, 
Wde < — 1 and Wdm > — 1; case II, Wde > — 1 and Wdm < — Ij since we should require that the final 
densities of dark matter and dark energy are both positive; as a matter of fact, the negative energy 
density may appear at most as a transient phenomenon, and it can not be a physically stable and 
permanent state. 

To study the stability property of the fixed points[21^, the evolution equations should be per- 
turbed to the first order and we list them in the Appendix. 

Then, the perturbation equations to the trivial critical point of [16] are simplified as 

2 

-^{SpdrnY = -Spdmi^ + ^^dm) + 5pdei-l + Wde), (19) 

2 

^{SpdeY = '5pdm(l - Wdm) " ^Pdci'^ + ^Wde)- (20) 

And near this fixed point, the eigenvalues of the linearized system read 

h = ^(-4 - 3Wde - Wdm - \J-8Wde + 9wj^ + SWdm " ^^WdeWdm + W^^), (21) 



h = ^(-4 - 'iWde - Wdm + \J-8Wde + 9wj^ + 8Wdm " WWdeWdm + wj^). (22) 

Stability requires that all of real parts of the eigenvalues are less than zero. Combined this 
requirement with the observational fact of Wde < we find 

- 1 < Wdm < 1, - TT-- < Wde < Wdm- (23) 

^ + Wdm 

A typical picture of the evolution is displayed in Figure 1. For illustration we set Wde = —0.4, 
Wdm = —0.1 and a = 0.9. As is seen to us, the orbits with initial values of energy densities 
lying in the first quadrant evolve to the second quadrant with the energy density of dark matter 



FIG. 1: The plane v versus u. The initial conditions are taken for different orbits. There is a stationary 
node at the origin of coordinates, which attracts orbits in the first quadrant; however the orbits have to get 
through the v-axis before reaching the node, and therefore the case may be unrealistic. 

negative, and they will keep negative before the tracks reach the trivial critical point at the origin of 
coordinates. Evidently, this case is physically problematic and might be irrelevant to our universe. 

Now let's turn to the non-trivial fixed point. Analytic eigenvalues of the corresponding linearized 
system to this critical point have not been found, and we will explore the stability property of the 
scaling solution via numerical method. Here we plot Figure 2 to display the properties of evolution 
of the universe controlled by the dynamical system. In the figure, we set A = 0, Wde = —1-2, 
Wdm = —0.2 and the positive coefficient a = 0.6 for illustration. The physically meaningful region 
is just the first quadrant and different initial value tracks evolve to the common non-trivial fixed 
point still in the first quadrant, which indicates the existence of the physically stable attractor. 

The most significant parameter from the viewpoint of observations is the deceleration parameter 
q, which carries the total effects of cosmic fluids. Using ([5]), ([8]), and ^ we obtain the deceleration 
parameter in this model 

_ _a J_ _ 1 /Odm(l + 3Wdm) + Pde(l + ^Wde) - 2A 
aH"^ 2 pdm + /Ode + A 

For a numerical example, we take the present dark matter partition uq = 0.3 and the present dark 

energy partition vq = 0.7, which is favored by present observations [1]. 

Figure 3 illuminates the evolution of deceleration parameter. As a simple example we just set 

Wdm = —0.2, Wde = -1-2 and a = 1.0. The dark matter dominates in the past, and since its 

effective EOS in that period is larger than the assumed EOS, dark matter decays into dark energy. 

Only some time ago, this decay made the dark matter stiffer and with the accumulation of the 

dark energy, the universe endures a deceleration-acceleration transition. The acceleration goes 



(a) 



(b) 



FIG. 2: The plane v versus u. (a) We consider the evolution of the universe. The initial conditions are taken 
for different orbits. It's clear that there is a stationary node, which attracts orbits in the first quadrant. 
(b)Orbits distributions around the node. 



— ^-r^rrTT^l ^ ^ ^ 0.5 

\ \ 
\ 1 

1 1 1 1 1 1 1 1 1 \ii 




-4 


-2 


\l 
\l 

-o\ 


\ 
.\ 


2 


4 



FIG. 3: The evolutions of q in the model (solid curve) and in ACDM (dashed curve), respectively. 



through the present era and, if the dimensionless coefficient a is about 1 or less, it will march for a 
phantom phase in the near future. That is to say, the deceleration parameter can be smaller than 
— 1. However, the super-acceleration will not last for ever. The decay process of dark matter into 
dark energy then reverses and dark energy begins to transfer its energy to dark matter component. 
This conversion impedes the super-acceleration and eventually the evolution of the universe is 



pulled back to the Einstein-de Sitter track; therefore, the troublesome cosmic big-rip is evaded 
From the figure one sees that current q ^ —1 and at the high redshift region it goes to 0.5, which 
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is consistent with current observations 



231 ] ■ As a comparison, the evolution of the deceleration 
parameter in a spatially flat ACDM is also plotted, in which the density parameter of dark matter 
^dm = 0.3 too. 

The density evolution of the cosmic fluid does not depend on the assumed EOS above, but the 
effective EOS. We define the effective EOS as the following procedure. Supposing the dark matter 
evolves adiabatically itself, we obtain its evolution from ([T]), 

da 

dpdm + ^Pdm+Peff) = 0, (25) 

a 

where Pe// denotes the effective pressure of dark matter. Then we obtain 

APeff ^[l-{l-2ay)wde] + 2^ + 3 + 2ay + {l + 2ay)wdm , 

Wdme = = ^— 1, (2bj 

Pdm 2 + lay 
^il-Wdm)-{l + 2ay + 3wde) + 2^ 
""'^^ = 2T^y ^' 

which is variable in the evolution history of the universe. Note that we have assumed both Wdm and 
Wde are constant from the beginning. Figure 4 displays the effective EOS of dark matter in which 
the same parameters are set as in Figure 3. It shows that, although its EOS is clearly negative, the 
dark matter behaves like cold dark matter until now. Only recently the dark matter has become 
stiff which relates to the current cosmic acceleration era. In Figure 5, different evaluations of 
the parameter a are taken and the corresponding deceleration evolution illustrated. A smaller a 
denotes a later deceleration-acceleration transition and a greater degree of acceleration in the near 
future. As the coefficient gets bigger, the cosmic evolution becomes closer to the concordant ACDM 
model and the transition happens more gently. For all that, the universe will go back to the same 
track of an everlasting de-Sitter phase with each of the energy densities unchanged thereafter. 



The deceleration-acceleration transition happens abruptly in the original semi-holographic 
model, and the present model in grand canonical ensemble is prone to moderate the process. 
As the chemical potential coefficient changes, the deceleration parameter tracks fill in the inter- 
mediate region between the original semi-holographic model and the ACDM model. Therefore, in 
this sense we can say that the particle decay picture extends the original semi-holographic model 
moderately and reconstructs the concordance model with the semi-holographic idea. The tracks 
may have or may not have a super-acceleration phase which depends on the different values of a. 
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FIG. 5: The evolutions of q in semi- holographic model with different values of chemical potential coefficient 
a [a = O.I, 0.6, 1.2 from thin to thick curves)and in ACDM (dashed curve), respectively. 

IV. CONCLUSION AND DISCUSSION 

Based on the semi-holographic model inspired by holographic principle, especially the previous 
studies of the relation between thermal dynamics and general relativity, we find that, the semi- 
holographic dark energy model in the grand canonical ensemble identification recovers the standard 
expansion history and evolves to the Einstein-de Sitter final state through a possible transient 
super-acceleration. 

A phantom phase may appear in the near future without a big-rip in the end. Assuming 
the total matter in a comoving volume to be evolved adiabatically, a compensating dark matter 
component should exist and the non-gravitational coupling between dark energy and dark matter 
is also compulsory. As is shown in the above, this interaction yields a future attractor solution, 
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which is a stable scahng solution for the dark matter-dark energy system in proper region of the 
parameters {wdm, Wde)- The final ratio of dark matter to dark energy only depends on Wdm, w^e, 
and is independent of the initial values of the densities of dark matter and dark energy and the 
coefficient a of chemical potential. This result is helpful to address the coincidence problem. The 
conversion of canonical ensemble context to grand canonical ensemble one seems to deform the 
original model to the well established ACDM model in the sense that different tracks roughly fill 
in the zone between the two models in the q — In a plane. Therefore, the particle decay picture 
extends the original semi-holographic model moderately and reconstructs the concordance model 
with the semi-holographic idea. This grand canonical ensemble viewpoint opens the possibility 
to deepen our understandings of dark energy thermodynamics and statistical properties, and is 
closely related to previous systematic application of the cosmological thermodynamics. 

With numerical examples, it is illustrated that the deceleration parameter is well consistent with 
observations. As a phenomenological model, the parameters including a may be further constrained 
by observational data and detailed investigations on global fitting of cosmological parameters will 
also be interesting. 
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V. THE APPENDIX: THE PERTURBATION EQUATIONS 

The perturbation equations for the combination {pdejPdm) are: 



de 



Z^Pdm 



-^Spde 



^ , Pde{l + 2aY + ?,Wde)^-pdm{l-Wdm)^ 



l + aY 



+ 



(1 + oYf 



PdeY + 1 + 2Q:1^ + "^Wde Pdm{'^ - Wdm)^ " Pde(l + "^CnY + 3Wde)a 



1 + aY 



+ 



2Y{l + aYf 



,(28) 
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3 + 2aY + (1 + 2aY)wdm + Pdm(l + Wdm)Y + PdeWdey 

1 + dY 



Pde[l - (1 - 2aY)wde]4. + Pdm[3 + 2ay + (1 + 2ay)^^rf, 



a 

mi 2Y 



(i + dy)2 

PdeWdey + /Odm(l + It'dm)^ + 1 - (1 - 2aY)'Wde 



1 + ay 



/)de[l - (1 - 2aY)wde]S + Pdm[3 + 2aY + (1 + 2ay)^i;rf, 



mj 2Y 



(l + aYy 



(29) 



with all the notations declared in the main text. 
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